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Abstract. We consider a one-dimensional optical lattice of three-dimensional 
Harmonic Oscillators which are loaded with neutral fermionic atoms trapped into two 
hyperfine states. By means of a standard variational coherent-state procedure, we 
derive an effective Hamiltonian for this quantum model and the hamiltonian equations 
describing its evolution. To this end, we identify the algebra C of two-fcrmion operators 
-describing the relevant microscopic quantum processes of our model- whereby the 
natural choice for the trial state appears to be a so(2r) coherent state. The coherent- 
state parameters, playing the role of dynamical variables for the effective Hamiltonian, 
are shown to identify with the ^-operator expectation values thus providing a clear 
physical interpretation of this algebraic mean-field picture. 
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1. Introduction 

Recently, a hierarchy of Hubbard-like Hamiltonians has been proposed to describe the 
behavior of ultracold fermions in one-dimensional optical lattices 



These lattices can be realized with a pair of lasers propagating at a given angle 
9 (8 = 7r represents the familiar counterpropagating case), with global confinement 
ensured by a magnetic trap (see Fig. [H for a detailed description of this setup see 
(J2)). The pair of lasers give rise to a interference pattern needed to obtain a periodic 
potential by AC Stark effect. The lattice constant can be adjusted tuning the angle 6 
according to the relation d = A/2sin(#/2) where A is the laser wavelength. In addition, 
it is possible to control both the barrier height of the periodic potential (as a function 
of the laser intensity) and the interaction between fermions via an external magnetic 
field (Feshbach resonance, see, e.g. (0)). 



These simple considerations allow one to understand how ultracold-atoms physics 
offers the possibility to explore experimentally a wide range of parameters set that 
would be unattainable in other contexts, such as the Hubbard model in condensed 
matter physics. 

As a first step towards the description of the cited models, we propose here a 
mean-field algebraic approach based on coherent-states procedure (0) for a fermionic 
one-dimensional array of harmonic wells. Although the analytical approach followed here 
may be regarded as completely general, future numerical analysis will concentrate on a 
dimer with a six-level structure per well as depicted in Fig. |2J). The approach followed 
here allows one a straightforward reformulation of the the usual mean-field approach for 
quantum system (based on the 'linearization' of the Hamiltonian and the subsequent 
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Figure 1. Sketch of the experimental setup considered. Two laser beams, propagating 
at an angle 9, give rise to an AC Stark induced periodic potential. In addition the 
fermionic cloud is confined by a cigar-shaped magnetic trap with trapping frequencies 
fl± and Q x . 
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solving appropriate self-consistency equations) in terms of a corresponding classical 
dynamical system. While, for fermions, the interpretation of the aforementioned 
classical dynamical system as a semiclassical approximation seems not beyond need 
of justification, it is possible to give a precise physical interpretation to the dynamical 
variables of the classical problem in terms of expectation values of quantum operators. 

In general, it is possible to consider a mean-field approach to a given problem as the 
constrained minimization of the Hamiltonian H over a algebra C. A different choice of 
£ will lead to different mean-field solutions (gQQ|). In particular, we will focus on the 
so(2r) coherent states that, as it will be shown, will lead to the Hartree-Fock-Bogoliubov 
(0) mean-field approximation, whose effectiveness has been proven for a single spherical 
harmonic trap in (Q). 

The paper is organized as follows. In section[2]a brief discussion of the general model 
considered will be given, along with some possible approximations in different physical 
situations. As we already mentioned, the fully-analytical control over the physical 
parameters allows to conceive various Hamiltonians that may have direct experimental 
relevance. In section El so (2r) coherent states and the relevant algebra will be defined. 
The end of this rather technical section will be devoted to the physical interpretation 
of the choice of so(2r) as the algebra for the mean-field procedure. In section 0] the 
classical Hamiltonian 7i c i will be deduced and the functional dependence in terms of 
quantum operators expectation values will be investigated. Finally, in section 14.11 the 
analysis of the classical dynamical system whose Hamiltonian is 7i c i is performed: Lie- 
Poisson brackets (namely the 'classical' commutators) and, consequently, the evolution 
equations for the dynamical variables are given. 

2. Model Description 

In (0), along the lines introduced in (j^), a generalized Hubbard Hamiltonian has been 
introduced to describe the behavior of alkali-metal fermionic atoms in a one-dimensional 
optical lattice of oblate three dimensional (2+ ID) Harmonic Oscillators (pancakes) 



In Eq. ((TJ), a must be considered as a multiple index a = {i a ,n a = 0, J a ,m a , a a } 
whose origin can be traced back to the space(+local)-modes approximation. In this 
picture n a , J a and m a are the local 2+1 D Harmonic Oscillator quantum numbers, i a 
is the site quantum number and a a is the spin quantum number. In the following we 
will confine our analysis to situation where radial modes only are involved in the system 
dynamics (i.e. we will "freeze" the axial quantum number n a to zero). The validity 
of this assumption is guaranteed as long as the radial trapping frequency Q± is much 
smaller than the axial trapping frequency u x , i.e. Q±/u x <C 1. In this case the tunneling 
coefficient assumes the following form: 




(1) 




(2) 
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if we allow nearest-neighbor hopping only 

Ti ai ip = T8i a+ ijp , (3) 

where T is a known function of the external parameters. Another assumption concerns 
the two-body interaction term U^p^g, which is treated within the pseudopotential 
approximation, leading to a delta-like spatial dependence, thus excluding nearest- 
neighbor two-body interaction terms. If the fermionic nature of the interacting particles 
is taken into account we have 

Finally, we give the expression for the one-particle energy term which is essentially given 
by the single-particle energy of the 3D Harmonic Oscillator 



A; 



fax ( n a + - 1 + HQ± (2J a + 1) - T ia>i 



(5) 



2 

where Ti a ^ a represents a "hopping correction" to the single particle energy term. 

For the case considered in Fig. |21 the selection rules imposed on the two-body 
interaction term select three possible values on U a ^,^,& that can be classified as: i) lowest- 
level/ lowest-level interaction terms, ii) lowest-level/highest-level interaction terms, iii) 
highest-level/highest-level interaction terms. 



(i„.i/„../„.m„.(T«} = {0.0.1/2.1/2.±l/2} 

c " - = 



c 



{i„.7i„..;„.m„.fT rl } - {0.0. 0.0. ±1/2} 



Figure 2. Schematic view of the single-particle energy levels. Up and down arrows 
represent possible spin values. We have explicitly indicated single-particle quantum 
numbers for the left-hand-side well. Analogously, it is possible to define quantum 
numbers for the right-hand-side well. 



3. Coherent States 

Since in our mean-field analysis we would like to keep trace of the (possible) atom 
pairing, the most appropriate coherent-states algebra seems to be, according to (4), the 
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algebra spanned by the r(2r — 1) operators {c£,cg(l < i,j < r), c a cp, c^ct}, i.e. so(2r). 
Its commutation relations can be written as 

\Xaf3i YyS] = |X* 2 /3> Y^ s ] = 

[ Y a/3> rfs] = Y^p5 a S + Y Sa Sfa - Y^ a 5p 5 - Y S/3 5, 

\Xaf3j Y^s] = Y^gSpy — YpgSa*, 

Y^ s ] = Yp 7 5aS — Y^Sps , 

having defined 



Y^p — C a Cp , 

Y 2 - r f r f 
Y^p — C a Cp . 



(6) 

(7) 
(8) 
(9) 



With the above definitions the coherent states can be expressed as 



10 > 



(10) 



CX]) - ^ (Va,f3clcl- H.C. 

To evaluate the expectation value of the Hamiltonian H defined by equation (|T|) 
over the coherent state of the form ([TUjl . it is necessary to evaluate the action of the 
operator 



Cl = exp - ^2 (vot^cicp - H.c. 

l<a^P<r 

over the fermionic raising and lowering operators. Namely 



C^c a Cl = exp 



~t *t * ~ ~ 



exp 



- (WH - rive]*) 



which, exploiting the BCH formula (jj) can be written as 



^ m 



E4-t _ * - - 
Vi,j C i C j Vi,j C j C i 



'J 



(12) 



(13) 



It can be shown that in the last summation the two first terms only survive, leading to 
to the following expression for Clc^Crf and CloyCtf respectively 

n f 4n = £ + 

i 

Clj^ C^Cl C-y ~\~ ^ ^Ta^^rn 
i 

with Qj = ZVij- 

We are now in the position to evaluate 

H d = (<f>\H-iiN\<l>), (14) 
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where the term fiN has been added to take into account the particle number constraint. 
With Eq. IfHty . Ea. lfI3jl becomes 



Hci = <0|O 



Ho + Hz 



(15) 



where 



Ho = 5^ r a>jB ct 



a,/3 

I 



#1 = U a,P,lA$f&°l 



with r Qi/ g = \ a 5 a ,(3 — Ta t /3 — fi5 a ,/3- Since Q is a unitary operator, we can write 



= u aJ 3 t vjttclpttcln£i%£i£i%£i. (16) 

The following expectation values must then be evaluated. For the one-body term 



a,f3 

and for the interaction term 



-a ^ ^ Ci, 



cp 



(17) 



C5 



C 7 + C; 7 c| 



4 + zJ 



(18) 



As it can be directly verified, in the calculation of the expectation values over the 
vacuum state |0) only the following terms survive 

(opH n\o) = EE r ^^(o|^io) , (19) 



a,/3 ij 



\i,j,k,l 



+ E^c 7 (oic4e4io)j 

i,3 / 



The two expectation values over the vacuum state give 
(0\c lC ]\0) = S i:i 

(0\ciCjclcl\0) = 5 u 5 jk - 5 jk 5u 
(0|cicjc 5 ct|0) = 5 if3 5 jS 

hence 



7~L d - E r «/3 E C* Q 0/3+ z2 Ua >P^> 5 

a/3 i a,P,l,5 



Y2i i^iaCjpCjsCi-y CiaCjpCi&Cj~i) + CpaC 



S'y 



i,3 



(20) 

(21) 
(22) 
(23) 

(24) 
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4. Effective Hamiltonian 

Hamiltonian (|2*4"j) can be shown to represent the effective Hamiltonian associated with 
H within the time-dependent variational principle procedure (0). The latter is based on 
approximating the quantum states of the system by a trial state 1 \P) satisfying the weak 
form of the Schrodinger equation (^/\ihd t — H\^f) = 0. Here, we assume that up to 
an irrelevant phase factor, is the coherent state defined in equation (fTU|) . The variational 
procedure allows one to derive the effective Lagrangian S = ($\ihdt — H depending 
on dynamical variables ( a/ 3, which in turn supplies the effective Hamiltonian (J2U). Such 
a procedure provides as well the dynamical equations pertaining to Hamiltonian ()24j) 
and the relevant Lie-Poisson brackets. The latter exhibit the same algebraic structure 
of commutators © and will be defined below. 

A quite direct physical insight about coherent-state parameters ( a/ 3 is achieved when 
considering the expectation values for the elements of the Lie algebra so(2r) over the 
coherent states We have 

(<P\c ] A\<P) = Q a , (25) 



\c%\ 



&C# = 6* > ( 26 ) 



(<P\c a cp\<p) = C a p, (27) 

showing how parameters ( a p are related to microscopic physical processes of 
creation/destruction of lattice fermions. Moreover Eq. ()24)1 can be written as 

7~Lcl = ^r Q/3 £ Q/ 3 + 2J U a fi }lj s [(ia-tips - €as£pr) + QaCs-y] ■ (28) 

The three terms represent the direct, the exchange and the pairing term in the HFB 
mean-field approximation. 

4-1. Evolution Equations for the canonical variables 

According to (J4J) the variables ( a/ 3, (^a an d £«/3 represent the canonical variable for the 



classical Hamiltonian H c i- With a well-known procedure (JlOj) it is possible to describe 
the time evolution of those canonical variables in terms of their Poisson brackets with 
7i c i. To write the Lie-Poisson brackets for the given dynamical system we can make 
explicit the structure constants for the so(2r) algebra 

C l,o,/3;l,7,<5 ~~ C 2,a,f3;2,-y,S ~~ U ' 

C l,a,j9;2,7,/' = ^<?,3 (8 m 8 v p8 a $ + 8 ^8^8 ^ — 8^8^8/38 — 8 ^8 v p8 al ) , 
C l,o,/3;3,7,/' m ' n = ^9,1 (^a8^s8f3y ~ 8^/38 u s8 a y) i 

C 2,a,/3;3,7,/' m ' n = ^9,2 (8^8 uy 8 a S — 8^ a 8„y8ps) ■ (29) 

Thus the Poisson brackets have the following form 
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{/> 9} — ^2 (^7/3^7 — 6ya^9<S + £<5a V7 ~ ^<5/3^a 7 ) ( 
a/376 \ 



9/ a/ 



+ 



9/ 




df dg 






d(,-yg d( a/3 






df dg 






dCjS d(* a/3 




dg 


df dg 









/> \ / df dg df dg \ 



+ (CaS&ya - £07^) ( ^ ^ - J • (30) 

Remembering that 

C P! e = {C p , e ,H} , (31) 
it is possible to write 

Cpfi = ^P a ^9a — ^OaCpa) + 

a 

/3 7 f? 
ceyrj 

/ y Cr? 7 {UpPyqipe ~ Ug plv ^j3 p ) (32) 
18717 

that provide the set of dynamical equations governing the evolution of the coherent 
state that approximates the system quantum state. In particular, they allow to find 
the mean-field ground the state for the system and to perform a weakly-excited state 
analysis. 



5. Conclusions 



In this paper we have formulated an HFB mean-field approximation for a one- 
dimensional array of oblate Harmonic Oscillators loaded with neutral fermionic atoms. 
As already pointed out by Grasso et al. (Q), the numerical solution to Eq. (|52~j) appears 
to be rather demanding from a computational point of view. It seems then appropriate 
for future work to concentrate on the simplest situation beyond known models like, as 
already mentioned, a dimer with a six-level local structure. 

In this case the evaluation of ground-state properties in this mean-field picture as a 
function of the relevant parameters(i.e T a> p, U a/ 3 7 s, /x) reduces to the fixed-point analysis 
of Eq. ()32j) . Moreover, an extension to finite-temperature properties does not seem 
beyond the possibilities of the analytical techniques here outlined and may represent 
one of the future lines of research. 



Mean-field algebraic approach 9 
Appendix A. Expectation values calulation 

In the present section we will explicitly calculate the terms obtained form Eqs. (jl7J18j) 
leading to Eqs. (fT9*j). To evaluate (fTTj) we need to perform the following product 
calculation 



af3ij 

which is equal to 



d p + Y 



a/3 



+ Y &&t> + Y ^ a c\ + ^ CLCjpCic) 



(a.i; 



(A.2) 



The evaluation of Eq. (JA.2j) over the vacuum state |0) leads to vanishing contributions 
for all non number-conseving terms and for all the terms with a lowering operator on 
the right-hand side (or a raising operator on the left-hand side) . Namely (see Eq. (fT9"|) 

<o|nt£ n|o) = EE r ^CO/3(o|Qc]|o). (A.3) 

a/3 ij 

With an analogous procedure it is possible to evaluate the expression given by Eq. 



G9 



a,/3,7,<5 



"a ~l~ ^ ^ Ci, 



7? + E C^ C J 



+ Y 



+ Y^4 



(A.4) 



leading to 



CaCpCfiC-y -\- 

Y c^44 a fcC 7 + ^44^ + 07^4444+ 

QpCaCjCS&y + E} Cj/3Cfc<54 a j4 a 7 + QpQ^C^CjCsc] + Cj/3 07^4^44 + 
Y C*a a i4 a 5 a 7 + E CaCHCi44 £ 7 + Y C*a0 7 a *4 a 4 + Ca^CfcSQC^cf + 
CiaCjfiCiCjCsC^ + (iaCjpCkSCiCjC-yCk + ^iaCj pClyC-i 

CjC S C\ + 



C*aC^Cfc507 a i a i a I a 3" 

With the same argument needed to obtain Eq. (|A.3|) we can write the expectation value 
of the operator defined by Eq. (|A.5|) over the vacuum state |0) as 



(A.5) 



(opH£\o) = U °>^ 



Y CaCK^o 7 (oi^4 a fio>+ 

\i,j,k,l 
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+ E^7<oi^ a 4i°> ( A - 6 ) 



h3 

which is the expression given by Eq. (f2T)|) . 
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